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o\ . 

We analyse the critical region of finite- (d)-dimensional Ising spin glass, in particular the limit 
of d closely above the lower critical dimension di. At criticality the thermally active degrees of 
freedom are surfaces (of width zero) enclosing clusters of spins that may reverse with respect to 
their environment. The surfaces are organised in finite interacting structures. These may be called 
■ protodroplets, since in the off-critical limit they reduce to the Fisher and Huse droplets. This picture 

provides an explanation for the phenomenon of critical chaos discovered earlier. It also implies that 
the spin-spin and energy-energy correlation functions are multifractal and we present scaling laws 
that describe them. Several of our results should be verifiable in Monte Carlo studies at finite 
0^ temperature in d = 3. 
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'. I. INTRODUCTION 

B 

'"^ ' No exact solution exists of the spin glass problem in the physical dimensions d = 2 and d = 3. 
C , The most successful approximate description was developed by McMillan 0, Bray and Moore 
Q ' and Fisher and Huse The authors speak of a scaling theory and the authors [||J|] of a 

droplet theory. In both cases the main hypothesis is the existence of an effective coupling constant 
dependent on the length scale L and satisfying scaling laws near zero temperature (with thermal 
exponent y) and at criticality (with exponent Uc)- These scaling laws may also be thought of as arising 
d ' from a renormalisation group (RG) for the spin glass, and in fact real space RG arguments are used 
by all the authors 0-^. In particular, RG arguments can be put forward |^ that experimental spin 
glass is well described by the Edwards- Anderson (EA) Ising spin glass model Q . 

Again, no exact RG is known for the finite-dimensional spin glass in dimensions d > 1. The most 
popular approximate RG is the Migdal-Kadanoff (MK) construction [0-0. In spite of its evident 
shortcomings the MKRG has been successful as an approximate calculation scheme for ferromagnets 
and other nonrandom systems. For the finite-dimensional spin glass there does not today exist a 
better alternative. 

Upon comparison one finds that the MKRG results for the spin glass are close, if not identical, to 
those obtained by the scaling and droplet theory. Most conclusions from MKRG are easily interpreted 
in the language of the droplet theory and the differences where they appear can be easily understood. 
There is therefore no reason at this level of approximation not to work with the MKRG, or, for that 
matter, with any real space RG of the same ilk. In this work we shall exploit the potential of the 
MKRG for spin glasses more fully than has been done hitherto. Our aim is to explain the phenomenon 
of "chaos" known to occur at criticality in low enough dimensions 12|. We achieve this goal by our 
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interpretation of new MKRG results. We are led, in particular, to a characterisation of the thermally 
active degrees of freedom in critical Ising spin glass just above their lower critical dimension. 

In order to set the stage for our work, we recall the most essential ones of the established results 
on zero temperature and critical chaos. 

A droplet of linear size L at a point R in space is defined as the domain of ^ L'^ spins around this 
point that can be reversed with respect to its environment at minimum free energy cost. The droplets 
are the low free energy excitations upon which the theory by Fisher and Huse is based. Droplet 
walls have an intrinsic width of the order of the correlation length and hence the droplet concept is 
naturally limited to length scales L > ^_(T). In particular, it no longer applies in the limit T Tc, 
where $-(T) diverges. 

Spin glasses have the property of being chaotic, as was revealed in the work of Bray and Moore 
1^, Fisher and Huse |^], and Banavar and Bray [ p^ . The term chaotic here denotes the property 
that the thermally averaged relative orientation of two spins a distance L apart is a rapidly and 
randomly varying function of temperature (see figure 1). If we denote the autocorrelation interval on 
the temperature axis of this function by AT^, then in the ordered phase 

AT^ L-'^ , asL^cx), (1.1) 



where C is the chaos exponent. Upon inverting (1.1) one finds 

Lat ~ AT~^ , as AT ^ 



(1.2) 



which is called the overlap length associated with a temperature difference AT. Numerically ^ takes 



values between 0.75 and 1 in 1 ^ c? 3 1 14- 16 r^,|l^ . Within the droplet theory the chaos phenomenon 
is elegantly explained as resulting from temperature changes that disturb the energy versus entropy 
balance of the droplet walls and thereby induce random droplet reversals. In this theory the sum 
2y + 2(^ of the zero temperature thermal and chaos exponents appears to equal the fractal surface 
dimension ds of the droplet walls. 

The original authors jl|-^ thought that chaos was limited to length scales L ^ ^_ in the spin glass 
ordered phase and to dimensions where such a phase exists, i.e., above the lower critical dimension 
dg. Ney-Nifle and Hilhorst ||ll|,|l^ subsequently showed that there is chaos also in the critical region, 
at least for spatial dimensions d in an interval extending upward from the lower critical dimension 
df w 2.5 P,P,p7t to some limit dimension estimated within MKRG to be d+ « 3.4 [ pl] , p^ . 

The reasoning can be summarised as follows. At a length scale L, the critical region is the 
temperature interval 

\T — T \ 

L—A ^ L-y^ (1.3) 



in which L ^ ^_(T). Here 1/yc = is the correlation length exponent. Ney-Nifle and Hilhorst use 
MKRG to calculate the effective couplings on scale L and find that inside the critical region they vary 
randomly with temperature with an autocorrelation interval 

AT^ - L-'^" , (1.4) 
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where the (positive) critical point chaos exponent is a new and independent exponent (see figure 1). 
This chaotic behaviour of the renormalised couphng (chaos with temperature) is reflected in the 
behaviour of the two-point correlation function < s^Si+L > as a function of temperature, which 
changes sign when the effective coupling at distance L changes sign. 



< SiSt+L > 



-1 -■ 



-1 -- 



This figure may be requested from the authors 



(a) L large 



+1 -■ 



-1 -- 



This figure may be requested from the authors 



(b) L still larger 



Figure 1. Temperature dependence of the two-point correlation function < SiSi+L > for (a) large L and (b) 
still larger L. At criticality AT^ ~ L'''" and near zero temperature AT^ ~ L^'' , for L oo. 



In view of this there is chaos inside the critical region if 

Cc>yc ■ 



(1.5) 



Figure 2 shows the dimension dependence of j/c and Cc- As d \ di, where de is characterised by 
y{df) = 0, these exponents join continuously their zero temperature equivalents y{d) and C{d)- The 
MKRG values for the critical exponents in dimension d — 3 are |ll^,0,n8| 



yc(3) « 0.36, Cc(3) « 0.57 



(1.6) 



Hence, according to these values, the 3c? spin glass is chaotic at criticality. 
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dg d+ 




Figure 2. The critical exponents yc{d) and Cc(d) as a function of the dimension d. They join continuously 
the zero temperature exponents y(d) and C,{d), respectively, as d approaches di from above. This dimension 
is defined by y{di) = 0. The intersection of the curves j/c(d) and C,c{d) defines d+. The curves shown are 
obtained from numerical evaluation of MKRG by the pool method O]. 



The problem that triggered the investigation of this paper can now be stated precisely. In view of 
the inapplicability of the droplet concept in the critical region, the natural question to ask is if the 
chaos at criticality can be given some different geometric interpretation. Blind analogy might lead 
one to suppose, for example, that the sum 2yc + 2Qc would be equal to the fractal surface dimension of 
objects yet to be identified whose reversal upon temperature changes would cause the critical chaos - 
however we shall see that this is not true. This work deals with answering the question raised above, 
but its results go beyond it: 

The paper is organised as follows. In section || we review and extend the existing MKRG theory 
of the Ising spin glass. Introducing the scaling of various averaged derivatives of the renormalised 
coupling constant , we obtain multifractal scaling of the energy-energy correlations at criticality. In 
section [II we show how the scaling laws deduced from MKRG can be given a geometric interpretation. 



The fundamental entities are surfaces (of width zero) enclosing clusters of spins that may reverse with 
respect to their environment. The surfaces are grouped into finite interacting sets that are self-similar 
in the range of scales between the lattice cut-off and their own size. Such a set of interacting surfaces is 
appropriately called a protodroplet , since in the off-critical limit it reduces to a Fisher and Huse droplet. 
Near the lower critical dimension there is confluency of zero temperature and critical phenomena. In 
order to clarify the picture, the characteristic lengths that play a role are analysed in section IV . In 
section ^ based on the physical picture of section III, we establish the multifractal scaling of the 
spin-spin correlation function at criticality. Furthermore, the effects of a magnetic field in the ordered 
phase and at criticality are discussed in a unified way, within the picture of section |l[ We use an 
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Imry-Ma type argument to deduce scaling laws for the leading behaviour of the magnetisation in the 
neighbourhood of the critical point, approaching Tc from below and above. The singular situation 
at di is also considered in detail. In section VI we point out the consequences of the strong spatial 
inhomogeneity of thermal fluctuations on the critical dynamics. The picture of this paper augmented 
with the additional assumption that the two-dimensional dynamics near zero temperature dynamics 
is activated, together with an argument of continuity, gives a first strong theoretical support for 
the conjecture that the critical dynamics of d-dimensional spin glass is activated in an interval of 
dimensions starting with di. In the conclusion, section VII, we mention possibilties of verifying our 
picture in simulations of the d= 2> EA model. 



II. EXTENDED SCALING LAWS FOR RENORMALISED COUPLINGS 



A. Introduction 



The Migdal-Kadanoff renormalisation transformation ||7[] was first formulated for the spin glass 
by Young et ai H,^. In this section we start from the same basic formulas. We briefly recall the 
main results known and then go on to extend some of them. Our interest concerns in particular 
correlations of various finds (with respect to the disorder average) between the renormalised and the 
unrenormalised couplings. 

The MKRG works as follows. For a rescaling factor h one combines a set of fe'' couplings at the nth 
level, {-fi"''}i„+i, to a new couphng, k\^'^^\ This is done by arranging the K^^'' of the block in+i in 
h''-^^ parallel strings (labeled by the index m) of b couplings each. [Whenever we will indicate a RG 
step explicitly in thwhat follows we will relabel the couplings of {-K^i"''}i„+i as -ft'^^.] The recursion 
relation then reads |^,|l9| 

K^"^ - E ^^^t^^l^ ri ^^^^^tl ■ (2.1) 

m=l fc=l 



In the remainder we set 6=2. The couplings of the originally given problem are of the form 
Ki = where i is a bond index, and the Ji are the coupling constants on the original lattice 

with root-mean-square (rms) value J. We shall also use K = Kf ^ = 13 J as the temperature variable. 
With the Ji given, the transformation (2.1) determines uniquely the couplings at all higher levels 
n = 1, 2, . . .. For the n times renormalised value of K we shall also write 

= when L = 2" (2.2) 



and for a specific scale-i coupling disc 

^« = ^1"^ (2.3) 



when in is the label of the coupling at R. K'^'^^ can be calculated from the unrenormalised rms 
coupling K = K'^^^ via the transformation (2J). If, as usual, after each RG step the distribution of 
couplings is replaced by a Gaussian with the same variance, one finds the recursion 
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arctanh^ J| tanhiiTj 



fe=i 



(n) 
fc 



(2.4) 



The transformation TZ has the foUowing properties (for some mathematically rigorous proofs see 
and references therein). It has the fixed points K — and K = oo, and, for dimensions d above 
di, also a critical fixed point if — Kc- Near the zero-T fixed point one has 



(2.5) 



with y < for d < di and y > for d > di, and near the critical fixed point K — Kc 

K^"+^^ - Kc-2y''{K^"^ - Kc) . (2.6) 



By expanding (|2.lD for /3 — > oo, one easily shows that closely above the lower critical dimension the 



equations ( |2.5| ) and ( 2^) can be written together in the compact form 



Co 



(2.7) 



where Cq is a positive constant. 



B. Multifractal scaling at criticality 



Let Kj^ be an individual renormalised coupling on scale L. One main result of reference [ pl] , p^ is 
that in the critical region K^ is a rapidly and randomly varying function of temperature with rms 
derivative 



/ dKL 
\ dK 



K-Kr 



Kc 



^L-y- 



(2., 



in an interval of dimensions de < d < d+. This chaotic behaviour of the renormalised coupling 
K^ is refiected in the behaviour of the two-point correlation function < SiS-i-^-L > as a function of 
temperature, which changes sign when the effective coupling at distance L changes sign (see figure 1). 



From (2.6) we have 



Upon rewriting (2_^) as 



(2.9) 



r dKli 

K^ ^ ~ L^-^ 

dK 



(2.10) 



one sees that K^ and ,r/^ are correlated, while for reasons of symmetry 
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(2.11) 



The relations (2.8) and (2.10) imply that, at criticality has a large random part, proportional to 



L^" , which has a random sign, and a small random part, proportional to L^'=, whose sign is positively 
correlated with that of (and which is therefore, in this sense, 'systematic'). 



Below we shall study the joint probability distribution at criticality of K'^ and ^Mf^Ki. We will 



argue that at criticality C,c is only the first of a family of an infinite number of exponents describing, 
together with ?/c, the scaling of the moments of the joint distribution. This is in contrast to the 
behaviour at zero temperature whose scaling is entirely described by the exponents y and C,. 



Since we wish to consider two-point correlation functions in section 11 C , we investigate the 
dependence of the renormalised coupling K]^ at a position R and at length scale L on the set of 
unrenormalised couplings {Ki] in the volume of order L'^ around R. The temperature derivative of 
can be calculated via (2.1) as 



dK 



E 



dKi K 



(2.12) 



where the sum is on the L"^ unrenormalised couplings. Neglecting cross terms (that contribute only 



to order L^^^ as -L ^ oo), we find from (2.8) and (2.12) at criticality 



-K, 



(2.13) 



Similarly ( |2.10| ) and ( |2.12| ) give 



(2.14) 



dK 

dK, 



The relations ( ^.13| ) and ( 2.14 ) are only two moments of a joint probability distribution of K^ and 



^Ki. For the other moments at the critical point we make the ansatz 



oK, 



(2.15) 



and 



(2.16) 



where ^c^'^ = 

Let no\ 
couplings 



Let now K^ = -f^f"^ and Ki = K'f^ and write by chain rule differentiation down the hierarchy of 



dK^^ " dK 



n 



(2.17) 
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where zi, i2, . . . denote the locations of the couphngs on the successive levels = 1, 2, ... to which Ki 
contributes. At criticality, the factors in the product are identically distributed random variables. If 
we assume them independent, then we have from ( 2.15| ) and (2.17) that at criticality 



K 



(9) 



(log 2) -1 log 



2q 



(2.18) 



The right-hand-side of this expression can be calculated from a single renormalisation iteration after 



which (^c'^^ is known ||20[]. Obviously there is no reason for (^c'^' to be linear in q. Hence we have found 
that at criticality there is multifractal scaling p^ . 

The above results are valid at criticality. It is instructive to compare them to their analogues in the 
ordered phase. There we have for L ^ ^_ [|ll|,|l2|J 



(9) 



dK 



L2y 



(2.19) 



where dg = 2( + 2y, hence 



{ OK, 



K, 



R dKi 



l2y-d 



(2.20) 



[We shall write dg throughout but note that this exponent takes the value dg = d—1 within MKRG.] 
This scaling is due to the fact that bonds contribute to the energy fluctuations if and only if they are 
cut by a droplet wall. In the limit T ^ 0, the derivative equals ±1 with probability iL'*""'' and 
vanishes otherwise [cf. (2.1) in the /3 cx^-limit, equation (3.3)]. For the same reason, we obtain for 
the higher moments: 



V 9Ki 



~dKi 



2q-l 



ld,-d 



(2.21) 



There is therefore no multifractal but only trivial scaling at length scales L 3> in the ordered 



phase. With the aid of ( 2.15 ) and ( 2.16 ) we can write the moments ( 2.21 ) in the scaling form 



2q 



2q+l 



(2.22) 



V dK, 

with scaling functions that satisfy 

•^I(O) 

^q{0) 
^qi^) 



X — >00 -7 



-2C 



(?) 



(2.23) 



a;2«-2/c-2C<'''+2Cc 
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In the limit d \ di along the critical line Tc{d) wc know that ?/c \ and Cc \ Ci^-i)- Since also ( 
2.15) and (2.16) must then reduce to (2.21), it follows that 



(2.24) 



It appears therefore that the zero- temperature exponent Cidg) fans out along the critical line into the 
family of the d'^K 



C. Energy fluctuations 



We shall now exploit our extended knowledge of the MKRG to study energy-energy fluctuations. 
These are, generally, obtained as the appropriate derivatives of the logarithm of the partition function 



logZ = ^G( 



(2.25) 



i/=0 



with respect to the couplings. Here expG*^'^ is the multiplicative constant that appears in the 
partition function in the vth renormalisation step. It is the contribution to the free energy PF coming 



from excitations on a spatial scale £ = 2"^ . Within the MKRG scheme defined by (2.1) the expression 
for G^") is 



(2.26) 



with 



m—1 



(2.27) 



where the sum in ( 2.26| ) runs over all blocks in+i of vi\\ level couplings and we have used the convention 
introduced in section ^ to relabel the in an explicit RG step. 



Let 



Ae^ = e; - < e., > 



(2.28) 



where Ci — SiSi', and Si and s^' are neighbouring spins in the original lattice linked by a coupling Ki 
The simplest energy-energy correlation function is then 



/32r2(ry) = < Ae,Aej >= ^ ^^^^ 



(2.29) 



A fourth order correlation function of interest is 



p'Uin,) ^ (3'[< Ae^Ae? > - < Ae? >< Ae| >] = 



i/=0 



(2.30) 
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These correlation functions are random quantities. We are now faced with the somewhat technical 
task of calculating averages as well as the averages of their powers. The resulting expressions are 
( ^.34 ), (2.37), and ( 2.38), and their derivation is as follows. Let for simplicity i = io and ii,i2, . . ■ 
denote the locations of the couplings on the successive levels i/ = 1, 2, ... to which Ki contributes and 
similarly for Kj . Typically, the couplings Ki^ and Kj^ fall in the same block when i/ attains a value 
v*{i,j) such that 



2" 



(2.31) 



where is the distance between the two spins concerned. They then contribute to the same 



and to the same K 



^ . Working out ( 2.29), one finds that the second derivative leads to 
terms of various types. The guiding principle in finding the dominant contribution for large rij is that, 



according to ( 2.15 ) and ( 2.16 ), averages with even powers of derivatives are larger than averages with 
odd powers, while some derivatives vanish due to symmetry properties of g{{K[^^}i^j^-^) and k[^^^\ 
Therefore we get 



/3'r2(r,,) 



OK, 



dKi 



n 



dKi ,dKj 



^» dK, ^3 



(2.32) 



Let /3^r2 (rij) be the term with v = n in ( 2.29| ), i.e., it is the contribution of scale L fluctuations to 



/? r2(ry ). Obviously, r2 (r) vanishes for L < r. Using again (2.16) we arrive at 



^ 2 



«-(i) 



1 



f.2d-2y^ 



(2.33) 



for the contribution of the energy fluctuations due to two spins a distance r apart in a volume of 
linear size L ^ r. Hence, summing over length scales, we get the energy-energy correlation function 
at distance r: 



1 



j,2d-2ya ' 



(2.34) 



the sum being dominated by the first term. 

If one works out the expression ( ^.30 ) for the fourth cumulant, the fourth derivative leads again to 
terms of various types. In order to find the dominant contribution for large we proceed as above. 
The result is 

P 141^1^ J _ oTK—gTK~^i \ dK, J \ OK, j 



K 



2 



3 \ dKj 



(2.35) 



The quantity [3 T^{rij), defined as the term with v ~ n in equation (p.3C|), vanishes for L < rij. For 



L ^r, averaging on the disorder and using (2.15) gives 



d-2Q 



1 



(2.36) 
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Integrating this equation over length scales L as above, one finds 



r4(r) = ^rf (r) 



(2.37) 



because d ~ 2C,c is the (necessarily positive) exponent of a moment of a probability distribution disc 
[cf. also (lITsI ) and ( ^^\. 

By taking powers of ( ^.32 ) and averaging with the help of ( 2.15 ) and ( 2.16| ), one finds 



(r) 



„4C<«' -2d 



(2.38) 



for the moments of the correlation fimction r2(7'). Analogous expressions hold for the moments of the 
higher order correlation functions. Since Cc'^\ given by (2.18), is a nonlinear function of g, this implies 



multifractal scaling of the energy- energy correlations at criticality. 

By summing finally ( ^.34| ) on all i and j in a volume V, one obtains the mean square fiuctuation of 
the total energy in that volume, 



< AE^ > 



Vd 



V 



drr''-'T2{r) ^ k^T^Cv{T,) 



V 



oo , 



(2.39) 



valid if d~2Qc > 0, which is the case in an interval upward from di [cf. ( p. 15 ) and (2.24)]. By summing 
(2.37) one gets similarly the fourth cumulant of the energy in a volume V , 



< AE'i > - 3< AE^ >2 



V V 

The result depends on the sign of 4^c — d and reads 



<Ag-*>^ 



V 



J drr'^-^r4ir) 



V^oo 




if ACc/d > 1 . 
if ACc/d < 1 . 



(2.40) 



(a) 
(b) 



(2.41) 



Equation (2.41a) is likely to hold in a narrow range of dimensions upward from di, and equation 



( 2.41 b) for higher dimensions. Equations ( 2.3E ) and ( 2.41 b) have nonanalytic correction terms in V. 
In nonrandom thermodynamic systems the energy cumulants tend to constants in the thermodynamic 



limit. Equation (2.41a), in contrast, expresses an anomalous scaling of the fourth cumulant with 
volume. It is easy to check that the same is true for the higher order cumulants. Such an anomalous 
scaling of the energy cumulants in spin glass was first pointed out by FH for the ordered phase [see 
[^, equations (7.11) and (7.12)]. The basic reason for this anomalous behaviour of the cumulants is 
that the energy fiuctuations are collective, even beyond the scale up to scale infinity, and involve 
larger and larger energy differences. 
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III. GEOMETRIC INTERPRETATION OF THE CRITICAL SCALING LAWS 



Our purpose is now to understand the scaling behaviour found in the previous section in terms of 
a geometric picture. In the ordered phase, for length scales L 3> the work of Fisher and Huse 
has made it clear how to pass from MKRG to such a picture. After recalling how this is done, in 



section Md, we will address the same issue in the critical region, in section [II B 



A. MKRG and the droplet picture of the ordered phase 



1. Droplets 



At length scales L ^ ^_ the block spins of the MKRG procedure are naturally identified with the 
droplets of the FH picture. Thermally active block spins, i.e., block spins coupled to their environment 
with energies less than k^T, then represent thermally active droplets. The reversal energy of a scale-L 
block spin is equal to twice the sum of the couplings that link it to its environment (with signs 
equal to the values of the neighbouring block spins). We shall speak of the as the excitation free 
energies of collective spin reversals on scale I/, even though strictly speaking it is always the sum of a 
few that is involved. 

The probability distribution of the scales with L in the same way as the probability distribution 
Pl of the reversal free energies AF^ of the FH droplets, VIZ. as 

1 /A \ 

P^(AF^)-^^(^) , ^^oo. (3.1) 

Here T(T) is the temperature dependent stiffness coefficient of the ordered phase with T(0) ^ J and 



T(T) |l6|,t3|. Equation (3J^) implies in particular that the probability f^{T) for a scale L 

block spin to be thermally active is 

Therefore in the large-i limit the thermally active scale L block spins become infinitely dilute and 
may be considered as independent degrees of freedom. 

FH argue that the droplet wall has a fractal dimension ds with d — 1 < ds < d. At low enough 
temperature the droplet wall will be of zero width and cut through a well-defined set of couplings of 
the unrenormalised system. The associated reversal energy is the sum of the energies of each of the 
couplings cut. The number of couplings cut will grow typically as L"^". MKRG gives dg — d — 1. This 



is easily seen if one expands the elementary renormalisation step (2.1) in the /3 — > oo limit. It reads 



= E (^il^il) li^ill) (3.3) 



m— 1 



at T = 0. Upon differentiating (^^) with respect to the 2'^ couplings that enter into the expression 
one finds that, at zero temperature, the derivatives with respect to 2'^~^ of them are equal to ±1, 
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whereas the others vanish. After n renormalisation steps with 2" ~ L one sees that the scale-i block 
spin corresponds to a droplet with wall area ^ L'^~^. This special value for ds, i.e., ds = d — 1, 
must nevertheless be ascribed to the particular choice of the MKRG recursion relations and is not 
fundamental. Numerical work ||2l|-|2^ points to nontrivial values of ds and one can envisage an 
improved MKRG that would lead to ds 7^ d — 1 . The wall of a droplet at finite temperature has an 
intrinsic width of the order of the correlation length (see figure 3). But since the RG flow carries 
the system to the zero temperature fixed point, the droplet walls coarse grained to at least the scale 
^_ have the same fractal dimension ds at length scales L ^ f _ . 



This figure may be requested from the authors 



Figure 3. A scale-L droplet at finite temperature {L ^ ^_). The wall of a droplet at finite temperature 
has an intrinsic width of the order of the correlation length In the large-L limit, scale-L thermally active 
droplets get more and more dilute and therefore may be treated as noninteracting. By the same token, it 
consists of one and only one main connected domain of order L'* spins, whose wall width is negligible with 
respect to its own linear size. Therefore droplets are two-level systems, with "up" and "down" referring to 
the main domain. 

In the large-L limit, scale-L thermally active droplets get more and more dilute and therefore may be 
treated as noninteracting. By the same token, upon assuming that there are no correlations between 
the spatial positions of thermally active droplets of different scales, one may neglect interactions 
between scale-L droplets and droplets on scales just below (e.g., L/2 and L/4). As one goes down 
further in scale, however, droplet-droplet interactions will at some point come into play, as may be 
estimated as follows. On some arbitrary linear scale £ < L, the wall of a scale-L droplet will pass 
through [L/tY' blocks of size i of which typically a number f^{T){L/ lY' will contain thermally active 
droplets interacting with the wall. For L below a value that following FH |^ we shall call Al, this 
number will typically be larger than unity and the wall of the scale-L droplet is "interaction widened" 
up to the scale A^. An expression for Al can readily be found but of interest here is only the 
property 

lim ^ = . (3.4) 

L — >oo L 

It implies that for large L a droplet consists of one and only one main connected domain of order 



13 



L'^ spins, whose wall width is negligible with respect to its own linear size. Therefore droplets are 
two-level systems, with "up" and "down" referring to the main domain. 



2. Interaction between droplets 



How does the interaction between droplets of different scales appear in the MKRG formalism ? To 
answer that question we consider a string of two scale-L couplings k["^ and ^^2"' (contributing to a 
higher level coupling k'^^^^) with in the center a scale-L block spin Sq and at the ends two larger 
scale block spins Si and 52 (one of them is necessarily a scale-2L spin and the other one belongs to a 



scale ^ 4L). The derivative of (3.3) with respect to a scale-L coupling is bounded by 



< 



dK. 



(n+l) 



dK. 



in) 
mk 



< 1 



(3.5) 



The preceding discussion implies that when its absolute value is not close to the lower or upper limit 
values or 1, respectively, this is the signature of thermal activity. If this is true for the derivative 



with respect to ifj;"-* or K2"'' , then the wall between the droplets Si and 5*2 may cut either bond 
of the string. The two positions of the wall differ by the orientation of Sq. Hence the droplet So 
thermally decorates the wall between the droplets 5*1 and S2- The two wall positions will in general 



have different weights. Upon working out the derivatives of interest with the aid of (2.1) we find 



dK. 



(«+i) 
+1 



isinh(2iff)) 



ok: 



(n) 



cosh 



K 



(«) 



cosh ( K 



,(n) 



K 



1,2 



(3.6) 



The two derivatives are of comparable absolute value, and hence thermal activity is nonnegligible, 
only when 



ifl"^ ± K 



(") 



J 



(3.7) 



The ± sign corresponds to the two possible relative orientations of Si and 5*2. When { p.7\ ) is satisfied 
for one of the two signs, it will not be, in general, for the other sign. This is the "interaction" between 
droplets of different scales. This feature of the MKRG also appears in the FH droplet picture: a scale 
i droplet lies on the wall of a scale L droplet (L > £), has a wall section of linear size i in common 
with it, and is thermally active for one, but not for the other orientation of the large droplet. 



At low critical temperatures the derivative (3.6) has a probability distribution on the interval [—1,1] 
which strongly peaks at and ±1 and is small elsewhere. At T = 0, it is the sum of three delta peaks at 
and ±1. Since the ordered phase is controlled by the zero temperature fixed point, the phenomenon 
of interaction widening discussed here plays no role for the large scale properties in the spin glass 
ordered phase. It is important, however, for the understanding of the crossover from zero temperature 
to criticality, as we shall see below. 
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B. The low-dimensional critical region - protodroplets 



1. Surfaces 



In the critical region, that is, on length scales L below the correlation length and whether below or 
above T^, the FH droplets cannot be defined. Yet, in the MKRG picture thermally active block spins 
can be defined as before, viz. by the condition that their interaction energy with the environment be 
less than some constant times k^T . 

Let us consider for definiteness the critical line in the T-d plane. On this line the probability for 
a block spin to be thermally active is of the order of k^Tcj J = p > 0, whatever the spatial scale L. 
Therefore there must be interactions between these spins whose importance is characterised by the 



parameter p. Similarly, from the fact that along the critical line the derivatives (3.6) have a nontrivial 
probability distribution it follows that there are interactions between thermally active block spins of 
neighbouring scales {L and 2L). 

Block spins can no longer be identified with droplets and wc should look for a new interpretation. 
The key to this interpretation is to consider the limit p — > 0. For d very close to di the scale 
independent probability p for a block spin to be thermally active is very small. Also in this limit the 



probability law for the derivatives (li) approaches the trivial three-peak distribution, which means 



that the interaction between block spins of neighbouring scales disappears. 

Let us now look at the unrenormalised lattice (scale L = 1) for p extremely small. The thermally 
active spins will be extremely dilute so that for the moment we may ignore the interactions in the 
discussion. We could indicate the thermally active spins in a figure by surrounding them with a 
surface. For L not too large we now identify a thermally active scale-L block spin with a size-L 
surface that encloses a cluster of unrenormalised spins that may reverse collectively with respect to 
their environment due to thermal agitation. These larger surfaces may be irregularly shaped. The 
essential point is that the surfaces of the individual thermally active spins as well as of the thermally 
active clusters are of width zero, because not interacting with any other thermally active degree of 
freedom. These thermally active surfaces are systems with two energy levels. They may be considered 
as independent degrees of freedom under certain circumstances, e.g., in the limit p <C 1 and/or for 



not too large values of L, or for the calculation of the spin-spin correlation function (see section V A) 



2. Interaction between surfaces 



We now ask about the effect of interactions between the surfaces. Interactions occur due to surface 
sections that they have in common. In the small-p limit, the thermally active surfaces occupy random 
positions in space. Let us consider two scales L and £ < L. If d{Tc) is the surface fractal dimension, 
then a size-L surface will pass through [L / 1)'^^'^"^ regions of linear size £ and interact typically with a 
number p{L/lY^'^-'^ of size t thermally active surfaces. This number becomes ^ 1 for 
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Hence beyond this length scale the surface entropy plays a role, i.e., the thermally active surfaces are 
likely to interact with smaller size thermally active surfaces with which they then form what we shall 
call composite thermally active surfaces. 

The probability for a thermally active surface to interact with another one of size comparable to 
its own is small (of order p), but is also taken into account at this level of approximation. We should 
therefore ask if the interactions do not lead to a coupling, direct or indirect, between all thermally 
active degrees of freedom in the lattice. It is not difficult to estimate that for small enough p the 
interactions couple the thermally active degrees of freedom only into finite disjoint sets. We call such 
a set an isolated thermally active composite surface (see figure 4). An appropriate name which is much 



shorter is a protodroplet , as will appear in subsection IIIC 



This figure may be requested from the authors 



Figure 4. A scale-L protodroplet in the critical region. The composite thermally active surfaces are 
many-level systems of coupled degrees of freedom. Since the composing surfaces may be of comparable 
size, the protodroplet cannot be simply characterised by a single up/down variable. It is a hierarchical object, 
composed of interacting two-level systems - the surfaces - which is statistically self-similar up to the scale of 
its own size. The energy fluctuations are due to thermally active bonds between adjacent spins that are cut 
by the surfaces, as indicated in the figure. The multifractality of the energy correlations is a consequence of 
the nontrivial distribution of relative weights between the two possible orientations of each composing surface. 

The composite thermally active surfaces are many-level systems of coupled degrees of freedom. Since 
the composing surfaces may be of comparable size, the protodroplet cannot be simply characterised 
by a single up/down variable. It is a hierarchical object, composed of interacting two-level systems 
- the surfaces - which is statistically self-similar up to the scale of its own size. Its multifractality is 
a consequence of the nontrivial distribution of relative weights between the two possible orientations 
of each composing surface. If the above argumentation, starting from thermally active surfaces and 
leading to the identification of protodroplets, is supplemented with appropriate scaling hypotheses, it 
may be seen as an alternative - valid in the limit d \ c?£ - to the MKRG formalism in providing a 
basis for the theory of the low-dimensional critical state. 
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C. Summary 



On the critical line Ti.{d) in the temperature-dimension plane the spin glass system is statistically 
invariant under scale transformations, and has thermal fluctuations on all length scales. The 
equilibrium ensemble is not dominated by a single representative configuration or class of 
configurations, which at first sight renders a geometric analysis forbiddingly difficult. Such an analysis 
is nevertheless possible, for sufficiently low Tc(d), i.e., for d close to d(. The simplifying feature is 
that at low temperatures, and whether the system is critical or not, the fraction of couplings that 
effectively participate in the energy fluctuations goes to zero. The thermally active couplings therefore 
determine a set of surfaces of low spatial density. Each surface encloses a cluster of spins which may 
reverse with respect to a frozen local environment in the ordered phase as well as at criticality. 



What has been said in section IIIB holds dXT = Tc as well as in the whole critical region. Imagine 
now that for some T below Tc we increase the length scale L. The probability for the protodroplet 
to be composed of two comparable size surfaces will go down from its value p for L ^ to zero for 
L — > oo. For L ^ ^_ a single surface will dominate and the protodroplet has become a droplet. Hence 
the concept of protodroplets generalises the concept of droplets in dimensions near dp and permits a 
geometrical picture that also applies to the critical region. Random reversals of protodroplets are at 
the origin of chaos in the ordered phase as well as at criticality. 



IV. CHARACTERISTIC LENGTHS NEAR THE LOWER CRITICAL DIMENSION 

Near the lower critical dimension there is a confluency of zero-temperature and critical phenomena. 
In order to clarify the picture, we analyse in this subsection the characteristic lengths that play a role. 



in particular the correlation length ^_ and the overlap length Lat defined in ( 1.2) 



A. Correlation length near 

The correlation length above and below is determined by the RG transformation TZ{K), whose 



explicit form near di, equation ( 2.7), can also be written as 



j^2L ^ 2^ ^ Co (K^)-^ + ■ ■ ■ , (4.1) 

where cq is a positive constant. For further analysis it is convenient to pass to continuous L, which 
gives 

yK"^ -ciK"^)-^ + ■■■ (4.2) 



d\ogL 



with c a positive constant. At low temperatures we may neglect the dot terms in (4.2). Above 



di, the exponent y{d) is positive and equation (4.2) has the two fixed-point solutions = oo and 
~ \fcjy = Kc- The latter relation can be recast as 



y = c ■ (4.3) 
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Since y is linear in d~ di, this shows that Tc ~ \/d — di just above dt |^^. At low temperatures, the 
solution of ( |4.2|) for initial value = iiT = is 



(4.4) 



This shows that just above di the critical point exponent ijc is equal to 2y, as also found by 
Correlation lengths ^+ (T) above Tc and ^_ (T) below Tc may be defined by 



(4.5) 



where /+ is a fixed constant often set equal to unity and /_ a fixed constant less than 1. These are 
the usual definitions; in particular, at {K^~)~^ = /_ the renormalisation group crosses over from 



critical to zero temperature fixed point control. From (4.4) and (4.5) we find 

\UJ) 



1 - 



- 1 



T>T,, 
T<T,. 



(4.6) 



Some comments are in place. Using 
well-defined limit 



), we find that for d \ de the coefficient of ^+ has the 



f+J 



log 1 



TT log 



(4.7) 



Hence the definition (4.5) fixes ^+ near di up to a finite multiplicative constant. At di, we have 



(4.8) 



as also found by McMillan and Bray and Moore If, however, we try to take the limit d\ de 



with ^ fixed in the expression for ^_ , we find 



lim ^_ (T) = lim 



rp2 



0, 



< ^ 



(4.9) 



so that the limit value of ^_ at ^ fixed is or oo depending on the arbitrary constant /_ . Of 
course, in the T-d plane, the point (T = 0, d = di) lies on the intersection of two lines, with ^_ = 
and ^_ = oo so that ^_(T, d) has a mathematical singularity in (0, c?^). Nevertheless this feature is 



insignificant if one keeps d fixed. At fixed d, equation ( [4.5| ) remains sensible. Equation (4.9) only 
indicates that the range of length scales L at which one is uncertain whether being below or above ^_ 
becomes larger upon lowering the dimension towards de- 
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B. Chaos near the lower critical dimension 



Chaos in the critical region has first been shown to exist by a RG analysis in which the density of 
zeros of a scale-L coupling on the temperature axis was considered. In this subsection we will illucidate 
in more detail the critical chaos near the lower critical dimension and show how it is reflected in the 
change in activity of the thermal degrees of freedom. 



1. The density of zeros of a scale-L coupling on the T-axis 



It is of interest to consider the autocorrelation function (T, T + AT) on the temperature axis of 
a renormalised coupling at scale L. For this quantity, normalised to 1 at AT = 0, one expects the 
behaviour 

1 - C^(T,T + AT) ~ iArVi(T) , as AT ^ , (4.10) 

where Pl(T) is the average density of zeros of such a coupling. This quantity is, in principle, 
determined when a renormalisation group transformation is given. 

Below we will calculate pl(T) near the lower critical dimension. Within the approximation of the 
Migdal-Kadanoff renormalisation group, Ney-Nifle and Hilhorst p^ , p^ found an expression for pl(T) 
which slightly rewritten reads 

TVi(T) ^ K^Y. {-^) M^(^'^)A(if(^"+^))...A(if ("-!)) . (4.11) 

3=0 ^ ^ 



Here, as before, K — the functions A and /i are determined by the renormalisation group with 
X{K) = 2^^^'^ increasing from 2^^ to 2^ in the interval {Kc^ co) and ^i,{K) behaving as ^jl{K) ~ ^ for 
K ^ CO, where ^oo is a constant whose MKRG value is w 0.12 . 



Passing to a continuous length scale and using (O) for the K'^'^ gives 



Moo 

\og2 Kl 



d£ 



1 U^y 



exp 




(4.12) 



In the limit of interest we can expand z(K) around K = oo yo- powers of K ^ . From equation (5.5a) 
of il2l we can deduce that z is even in K so that to lowest nontrivial order 



(4.13) 



Using (4.13) and again (4.4), and setting r -'^ = ^ — 1, we obtain 



1 \ t 



1 + 



(4.14) 
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Changing variables to s = ^ , we get finally for the average density of zeros on the temperature axis 



with 



2y 



dss ^ «+2c (1 + s)" 



_ 7_ <" 



(4.16) 



Pl{T) 



d\de 



2C log 2 



fejT 



We are now in a position to consider various limits. From (4.15) and ( 4.16| ) we get 



2Cc log 2 



I T-Te 



^ - 1 





< 1 



(a) 
(b) 



(4.17) 



The expression ( [l.l7t 3) shows that for large enough L the density of zeros has a negative derivative at 
Tc, and must therefore have a maximum between T — and T = T^. For smaller L, the derivative at 
Tc is positive and the maximum of pl(T) must occur for T > Tc- The crossover value is 



~ exp ' 



'l0g(l + 7^)^ 



22/ 



(4.18) 



In the limit d \ de at L and T fixed, we find 

d\di I /loo 



Pl{T) 



2Clog2V J3 



(4.19) 



valid in the whole interval of temperatures ^ T ^ Tc. The ^/T behaviour at low temperatures was 
known for d > dg H,|l^,|l^ and for d < dg |ll|,|l^. The present calculation shows that this behaviour 
persists when d passes through dg. 



2. Confluency of critical and zero-temperature phenomena near dt 



A last comment on the confluency of critical and zero temperature properties within the geometric 
picture near dg is in place. Each sign reversal of a coupling on the temperature axis will typically 
be accompanied by a sign reversal of the correlation < sqSl > between two spins a distance L apart 



in the region around R. From ( 4.19 ) for the total number of sign reversals of < s^sl > between 
and Tr 



Tl{L) - /drpi(T) 




(4.20) 
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valid near d^. So, typically, as one passes from the ground state to the critical state there are no sign 
reversals of sq with respect to sl below a length scale L^^ given by 

Hence below this length scale, a protodroplet around a given point R at criticality coincides typically 
with the droplet at zero temperature in that region, even though an average over the whole system 
yields critical scaling and statistical invariance with length scale. 



V. SPIN-SPIN CORRELATIONS AND RESPONSE TO A MAGNETIC FIELD 



A. Spin-spin correlation functions 



We wish to consider here spin-spin correlations at criticality. Instead of using the MKRG formalism 
as we did for the energy-energy correlations, we shall now make freely use of the interpretation in 



terms of thermally active surfaces developed in section IIIB . For the following calculation we consider 
the thermally active surfaces as independent two-level systems, i.e., we neglect their interactions. 



Let < SiSj >c 



-< SiSj > — < Si >< Sj > he the connected correlation function between two spins 
a distance L = 2" apart. In this section i and j denote spin sites and not, as before, bonds. By 
a "surface" enclosing a given spin at a given scale we mean in what follows the one with the lowest 
excitation energy of all possible such surfaces; it need not, and in general will not, be thermally active. 
With respect to the thermally active surfaces introduced above these surfaces are exactly what the 
droplets are to the thermally active droplets. 

Let k be the index denoting the surface of scale 2*^ that encloses spin i but not spin j, and let 
similarly I be the index of the scale 2' surface that encloses spin j but not spin i. Furthermore let m 
be the index of the scale 2™ surfaces that enclose both spins. We shall denote the excitation energies 
of these surfaces by e^, e;, and e„i, respectively. 

The spin clusters enclosed by the surfaces introduced here can now be treated in the same way as 
the FH droplets in the ordered phase |^,^ . For the spin-spin correlation it is sufhcient to work in the 
approximation of independent surfaces; taking their interaction into account would only renormalise 
several of the constants but not lead to any new effects. 



Upon treating the surfaces as independent two-level systems one finds 



Ji-i „ n-i „ 



k=l 



1=1 



(5.1) 



where s^Sj is the relative orientation of Si and Sj in the configuration with none of the surfaces excited. 
From this expression we obtain the averages 



< SiSj >29 



Yl tanh^' 



Y[ tanh' 



2q 



(5.2) 
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At criticality, all the e's in the products are independently distributed random variables. Therefore 



the last factor in (5.2) equals unity and we find 



Obviously, we have in general Tyc'"* 7^ QVc^"^ = QVc- Hence, the rjc'''^ are another series of independent 
critical exponents, i.e., we have also multifractal scaling of the spin-spin correlation function at 
criticality. 

In , Sourlas considered the long-distance behaviour of correlation functions of critical spin glass 
using an effective field theory and the replica trick. For d > 6, he found after averaging over the 
disorder that the different replicas are not coupled by relevant or marginal operators in the effective 
replicated Hamiltonian. It follows as in that 77c'''' = for d > 6, i.e., the mean-field (Gaussian) 
behaviour. This is not the case in d = 6 ~ e dimensions. An e-expansion around the effective field 
theory in d = 6 leads to multifractality, i.e., r]f^ ^ <lVc:^\ as found in this work by a very different 
approach. 



B. Response to a magnetic field 



1. Introduction 



In this subsection, we will consider the magnetisation in a small uniform magnetic field H . Many 
of the results obtained in this section can be found scattered in previous work by several authors 
§,|j2|J|. McMillan §], Bray and Moore |,|2|, and Fisher and Huse |,| analysed the effects of a 
magnetic field on the ordered phase of a spin glass by applying a variant of the argument used by 
Iniry and Ma |2^ to determine the stability of Ising ferromagnets to random fields. They find that 
the zero-field ordered states are destroyed beyond a length scale £,h given by 

with T(r) ~ J^Z^ and q{T) ~ qoC-- Equation (|]J) holds as long as the field H is small enough so that 
^ £.^- Here T(T) and q{T) are the temperature-dependent stiffness coefficient and the Edwards- 
Anderson order parameter, respectively. The corresponding magnetisation m{H) is calculated to be 

The behaviour of m{H) for T = Tc can then be obtained on general scaling grounds [^,0-^. 

The picture developed in this paper now permits us to extend the Imry-Ma type argument to the 
low-T critical point. The crucial property common to both the droplets of the ordered phase and 
the composite clusters at criticality is that they may reverse with respect to a fixed background . The 
effects of the magnetic field, in the ordered phase and at criticality, can therefore be summarised as 
follows: 
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• The field ahgns small composite clusters (thermally active or not) with a certain probability and 
larger ones with a larger probability and 

• disrupts the fixed background by aligning practically all composite clusters from a certain length 
scale (~ ^h) on. 

This allows us to present the effects of a magnetic field on the system in the ordered phase and at 
criticality in a unified way. 



C. A magnetic field at criticality 



1. Temperatures below Tc 




We will first give the expressions for the correlation length in a field and for the magnetisation 
m{H) for T ^ Tc, valid irrespective of the relative size of the variables t and h. In order to do that, 
we have to know the typical reversal free energy AF^ and the typical magnetic energy E^^^ = H 
of composite clusters, of linear size £, where denotes their total magnetisation. We recall that at 
length scales £ ^ these composite clusters are the FH droplets. The reversal free energies AF^ are 

AF'^i ' ' ■ (5.6) 

^ e_ ^ £ . 

In the ordered phase, the mean square value of the magnetisation is q(T)£'^ for large droplets. At 
criticality, we have to take into account the renormalisation of the order parameter 

QiiTc) ~ qor^^ . (5.7) 
Hence, for general T, the expression is 

(5.8) 

^- < e. 

From simulations (see |^ and references therein) and experiment one can infer ijc ~ —0.3. In 
both cases rjc is obtained by hyperscaling from the values of the critical susceptibility exponents. 

Flipping a domain of scale i will change the system's magnetisation by an amount of order £^ , 
the magnetisation of the domain going from to —M^. Half of the composite clusters can lower 
their total free free energy by aligning with the field. Thus, the system breaks up into domains of 
typically the size of a correlation volume The correlation length S^h is obtained by comparison of 
the reversal free energy and the magnetic energy as the solution of 

Ai^«" . (5.9) 

The typical fraction of domains of linear size £ ^ £,h that will reverse to align with the field is 
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mag 

AF<^ 



(5.10) 



Finally, the total magnetisation M (H) in a volume (,fj is given by the sum over all the domains, up 
to scale ^ff, that align with the field: 



(5.11) 



Hence the magnetisation (per spin) m{H) is 



M (M^ 



(5.12) 



Introducing ( |5.6| ) and (5^) into (5_^), and solving for S^h, we get for the correlation length in a field 
near criticality: 



(5.13) 



where A = v{d + 2 — ?7c)/2 is the gap exponent and hyperscaling is tacitly understood. Numerical 
values for A can be obtained with the help of the relation A = /3 + 7 where 7 is the exponent 
describing the divergence of the nonlinear susceptibility. From Monte Carlo simulations and high- 
temperature series expansions one has estimated (3 = 0.6 ± 0.1 and 7 = 2.95 ± 0.25 (see ||2^] and 
references therein). Analyses of experimental data give the same range of values for (3 but favor a 



higher value of 7 = 4.5 ± 0.3 |30[|. The scaling (5.13) at criticality has been obtained by Fisher and 
Huse |3j previously [their equation (5.13)]. The scaling function J-{x) satisfies 



J^(0) = est. , 

!F{x) analytic for x — Q 



(5.14) 



The magnetisation is obtained from (5.6), (5.8), ( |5.9| ) and ( 5.12): 

m{H) -xH = with x ~ d~-2+r)^ J 

and the singular part of the magnetisation m^g verifies 



(5.15) 



2?! 



t 



d-2 + T]c 

with a scaling function A4{x) such that A^(0) is a positive constant, and 



(5.16) 



Mix) 
M{x) 



analytic for x = , 

±^{~xr-^^ 



(5.17) 
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The limit x ^ — cxd is in agreement with the result obtained by Bray and Moore |2j [their equation 
(46)]. The linear term xH exceeds the magnetisation m{H) by the singular contribution m,g(t,H), 
because, just as in the ordered phase j^, the low-field linear response of the domains with L ^ has 
saturated. 



In particular, we see from (5.16) and (5.17) that the exponent of H in m^g is 1 + 2 and 1 + 2-^ 



for a; — *■ and x — *■ — oo, respectively. Hence, the singular contribution is dominated by the linear 
response. However, both values of the exponent are smaller than 3, so that the singular term dominates 
over the term in the magnetisation and causes the zero-field nonlinear susceptibility to be infinite. 



Furthermore, we see from (5.16) that the third derivative of the magnetisation diverges, for H ~* 0, 
at the critical point as 



dH^ d-2 + T]c 




1+2^ 



. (5.18) 



For dimensions d close above dg the following inequality holds: 



y <l^d<(3^. (5.19) 



d — 2y A vy 
The second one follows from the fact that, for d\ d^, we have i^y ^ i Hi y \ and /3 ^ y~i /" oo 



(see (^). It seems reasonable to assume that (5.19) is valid in general dimension d > de, so that the 



divergence for iJ ^ of-|^ is stron ger in the ordered phase than at criticality. 



Summarising the above results for H ^ at t fixed, we get: 

T<Tc: m{H) = awH - as{t)H'^+^ + 



a,(t) ~a,o|i|^ ^ (5.20) 



T = Tc: m{H) = a^H - a''^H^+^ + ... 

Here, aio, a^o, and are amplitudes, and ( ^.19| ) implies that the exponent of as{t) is positive. 



2. Temperatures above Tc 



Let ^+ be the zero- field correlation length for T > Tc. Domains of linear size £ ^ ^+ behave 
critically, whereas much larger domains are paramagnetic. When H is small, domains of size up to 
^^ji align partially and in order to obtain the magnetisation we have to integrate until scale f + . When 



H is larger, there is < £,+ , and the integration runs until (^h as in equation (5.12). In the latter 



case, we get again the result (5.16), but with small positive argument of the scaling function. In the 
former case, we have 



€+ 

m{H) ^ H j 



d£ {M^ 
£ £dAF^ 



(5.21) 



so that 
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m{H) - xH 



with X 



d-2 + r]c J d-2 + ?]c J 

the singular part of the magnetisation contributing a term ~ to the hnear susceptibiUty. 



(5.22) 



We can write the magnetisation for T ^Tc, from ( 5.16| ) and ( 5.22), in a compact form: 

m{H) — xH • 



(5.23) 



with a constant contribution of the nonsingular part of the magnetisation to the hnear susceptibihty 



as in (5.22), whereas the singular part of the magnetisation behaves as 

m,g(t, h) 



d-2 + r]c \ht 
with A4 (0) a positive constant and the scaling function A4 {x) satisfies 



M{x) 
M{x) 



analytic for a; = 



(5.24) 



(5.25) 



To study the limit ft, ^ at small positive t of (5.25), we write first 



\h'S 



where M = (y^r) • From ( we have 

M{u) ~ , 



(5.26) 



(5.27) 



For t > 0, the function m^g should actually be nonsingular in h. Hence, with the assumption that M 
is analytic in h'^, we get in summary for the magnetisation above the critical point: 



T>Tc 



m{H) = ai{t)H + a3{t)H^ + 
ai{t) ~ aio - aipf^ 



(5.28) 



where aio, ai/3, and ai^ are amplitudes. This shows that, for t \ 0, the linear susceptibility x 
has a nondiverging singular contribution ^ and the well-known feature [ p5[ that the nonlinear 
susceptibility diverges as t~'^ . 



3. A magnetic field at d = di 



We now turn to the special case of applying a small magnetic field H at the lower critical dimension 
d — df . For d \, di, the contribution of the nonsingular part of the magnetisation to the linear 
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susceptibility [see (5.15)] diverges, because rjc 2 — d. However, due to the simultaneous divergence 
of the singular part of the magnetisation [see (5.16), ( p3.24 )] the linear susceptibility remains finite 
for r > 0, because the leading contributions of the nonsingular and of the singular part cancel each 
other in the limit d = di, as previously observed by [|l|,|25||. This gives rise to a logarithmic singularity 
of the magnetisation, for _ff at T = 0, which can be seen as follows. To obtain the expression 
of the magnetisation at d = di, we can reason as in the preceding subsection for T ^ Tc, with here 
Tc{dg) — Q. At T = 0, the reversal energy AF^ and the magnetisation of a cluster of linear size £ 
are 



J , 



(5.29) 



for all length scales £, so that from ( p.9\ ) 



J 



(5.30) 



However, for T > 0, the renormalisation of the reversal free energy AF is nontrivial. Solving (4.2) 
for y = we obtain 



/^F'^J^^l-c{^-^\ log^ 



(5.31) 



as McMillan did previously [|]. The magnetisation, for iJ ^ at T ^ 0, is readily calculated in the 



same way as above for d > d( and T ^Tc, the zero- field correlation length being ^+ ~ e ' 



see ( 4.8)]. One finds 



m{T,h) = ^^Md,{k^T\ogh-') , 



(5.32) 



with a scaling function Mdi that satisfies 



Md.ix) 1 



(5.33) 



Thus, in the limit H ^ at T fixed, we obtain the Curie law as expected in the paramagnetic phase 
on general theoretical grounds [pTI]. 



VI. SOME COMMENTS ON CRITICAL DYNAMICS 



Random systems, like random field magnets and spin glasses, have far slower dynamics than 
conventional pure systems (see e.g. p^-p^ and references therein). This is often attributed to fact that 
the competition between various types of interactions leads to free energy barriers of height that 
have to be crossed in thermally activated dynamical processes. If a static theory like the one of the 
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preceding sections is to be extended to include dynamical phenomena, then additional assumptions 
concerning these barriers have to be made. It is generally assumed that a barrier typically grows with 
the length scale L involved in the process as 

^L'l' , (6.1) 

with a harrier exponent > and gives naturally rise to long time scales 

- exp(^) . (6.2) 

This type of scaling is called activated dynamic scaling. The activated processes can be either 
fluctuations in equilibrium or involve approach to equilibrium. 

Generally, it has been agreed upon that this is indeed the mechanism valid at sufficiently large 
length scale in the spin glass ordered phase. However, it has remained unclear whether the critical 
dynamics in these systems should be described by activated or rather by conventional (power-law) 
dynamic scaling. Until now, most often critical dynamics has been scaled in a conventional manner, 
with anomalously large exponents, although thermally activated critical dynamic scaling with a small 
barrier exponent also appears to account reasonably well for experimental spin-glass behaviour at the 
critical point in three dimensions and does not seem to be incompatible with numerical simulations 
P^ , ^ . Again, however, the scaling with either approach can be corrupted in d = 3 due to a number 
of serious problems (see [ ^0|j29| ] and references therein for details). 

The picture of this paper represents an extension of the droplet picture of the spin glass into the 
critical region. We argue in the following that it suggests that spin glass dynamics is activated also 
at the critical point, at least in dimensions close above dg. We recall that in the FH theory the 
droplet walls minimise the free energy of a reversed spin domain. They are locally optimal and the 
fact that they dominate the physics of the ordered phase reflects the strong spatial inhomogeneity 
of thermal fluctuations in that phase. In particular, to reverse a droplet of scale L, the domain wall 
must pass through regions of high free-energy cost which are larger than the excitation free energy 
of the completed droplet. Therefore FH argue that droplet reversal is a thermally activated process 
with barriers B^ behaving as in (6.1) and growing at least as fast with length scale as the free energy 



difference between initial and reversed state. Therefore its barrier exponent satisfies 

y{d) < i,{d) (6.3) 



in dimensions d above di |3j. 

The same reasoning can be applied to the protodroplets of low-dimensional spin glass at criticality. 
But the lower bound ^ V-'c for the scaling of the barriers at criticality, obtained from comparison 
with the scaling of the free energy, is trivial. Indeed, -00 = would mean conventional scaling. 
We now give an argument that points towards ipcid) > in an interval of dimensions extending 
upward from di. Dynamics near zero temperature in two dimensions has been argued on the basis 
of numerical simulations and experiments to be thermally activated. Numerical work points to 
-0(2) w 1 (but see also [Q and references therein) and experiments give ■0(2)/|y(2)| — 1.6 ± 0.2 and 
ij{2)/\y{2)\ = 0.9 ± 0.1 for 40A ^ and < 20A Q thick spin glass films respectively. This would 
mean ipid) > 0, and, since the barrier exponent is expected to increase with dimension, also ip{di) > 0. 
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Hence by continuity ijjc{d) > in an interval of dimensions upward from . In these dimensions the 
barrier heights scale at criticality like 



B^{T,)r^ JL^' , 0<ij,. (6.4) 

Hence, the characteristic reversal time of a composite cluster at a large scale L at criticality is 
exponentially length scale dependent: 

~ To exp — — . (6.5) 

Such a behaviour implies j^] that ac correlation functions should scale, in the critical region, as 
functions of '"^^^"J"-* and that 

X3(3lj;lj)- \ \og^ {ujtq)\ , (6.6) 

i.e., the cubic nonlinear susceptibility diverges only logarithmically with vanishing frequency at 
criticality. 



VII. CONCLUSION 



We have analysed the critical state of the low-dimensional Ising spin glass, in particular in the limit 
of the dimension d approaching the lower critical dimensionality di. 

The best existing approximate theory for the low-dimensional spin glass is the scaling or droplet 
theory developed since the middle 1980's by McMillan, Bray and Moore, Fisher and Huse, and others 
The building blocks of this theory, in the language of Fisher and Huse, are the low free energy 
excitations, called droplets, from the spin glass ordered state. Droplets are reversed regions of one 
phase embedded in the opposite phase. Their excitation energy has a random value on an energy scale 
that increases with their size. The droplet wall, like in fcrromagnetism, has a width of the order of 
the correlation length. Consequently droplets exist only in the ordered phase and the droplet theory 
in its original formulation says nothing about criticality. 

We show in this work that the FH theory has a natural extension into the critical region, that is, to 
length scales less than the correlation length. This extension can be presented as a heuristic expansion 
around the dimension di. Its main idea is summarised as follows. 

In a critical state very close to di most bonds are frozen, but there is a very dilute set of surfaces (of 
width zero) that enclose clusters of spins able to thermally reverse with respect to their environment. 
The density p of these surfaces is the expansion parameter; it is proportional to the critical temperature 
Tc. Because of scale invariance there is the same density p of thermally active surfaces at all length 
scales. Since p is small but nonzero, interactions between these surfaces (through common surface 
sections) should be taken into account. Therefore, in the low- dimensional critical state the thermally 
active degrees oj freedom are organised in finite disjoint sets of interacting thermally active surfaces. 
These sets have a distribution of sizes that ranges to infinity. Each set has a self-similar structure on 
scales from its own size down to the lattice cutoff. Such a set of interacting surfaces is appropriately 
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called a protodroplet. In the off-critical limit, that is, for length scales beyond one of the surfaces 
participating in a protodroplet will dominate in size so that all the others may be considered as 
thermally decorating this large one and building up its wall: the protodroplet has then become a 
droplet. 



This picture of the critical state explains the phenomenon of critical chaos |]ll|,y_2|. It has enabled 
us to discuss the effects of a magnetic field in a unified way, using an Imry-Ma type argument in the 
ordered phase and at criticality. Finally, we have briefly argued, on the basis of additional hypotheses, 
that in low enough dimensions the critical dynamics is thermally activated. 

The picture of this paper has the practical advantage of making definite statements about spin 
glass properties on length scales below the correlation length (in contrast to the Fisher and Huse 
droplet picture, which is valid on scales above the correlation length). Some of our predictions should 
be easily accessible in simulations of the Edwards-Anderson model. One may of course attempt to 
determine the existence and properties of the thermally active surfaces (they exist from the scale of the 
lattice distance up) , but also more conventionally verify our predictions for various spatial correlation 
functions. 

The following three remarks about possible future sinulations are in no way meant to be exhaustive. 

1. About thermally active cluster surfaces. A quantity that measures the degree of thermal activity 
of the surface element perpendicular to a bond between two neighbouring spins Si and Si' is 
I < SiSii > \. If this quantity is zero, the bond is strongly thermally active; if it is close to 
unity the bond is nearly frozen. Its spatial (or disorder) average is equal to 1 — ©(fcsTc/J). 
Conceivably a chart of all nearest neighbour | < SiSi' > \ could be converted at the end of the 
simulation into a map of all thermally active surfaces. 

2. Dimension dependence. As the dimension d moves up from dg, the critical temperature Tc moves 
away from and thermally active cluster surfaces become gradually harder to distinguish, since 
almost all bonds will acquire some degree of thermal activity. Also, the thermally active surfaces, 
even if they can be distinguished, may, from some dimension on, form infinite percolating 
structures and the theory of our paper is not directly applicable. In dimension d — 2 and 
at sufficiently low temperature, the thermally active cluster surfaces exist on length scales less 
than the correlation length ^, but the correlation functions are not multifractal; beyond ^ the 
surfaces form an infinitely connected structure. Therefore d=3 is the ideal dimension for testing 
our predictions. 

3. Analysis in terms of correlation functions. This is the more conventional approach. The theory 
predicts power laws at distances less than the correlation length (but of course sufficiently large 
with respect to the lattice cutoff; how much larger we do not know). The powers are in general 
different for different moments of a given correlation function. We extract from the preceding 
sections what seem to us the most important ones. 



The disorder averaged energy-energy correlation function < Ae^Aej >, equation (2.28), should 



behave as (see equations ( 2.29) and (2.34)) 



< AciAej > ^ \i - jl^y^'^" {d>de, \i - j\ < . (7.1) 

However, the theory also says that < Ae^Aej > has fluctuations larger than its average value. 
Its variance (see equations ( 2.30| ) and ( 2.37| )) is 
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(The Migdal-Kadanoff values for d = 3 are « 0.36 and (c « 0.57 [|ll|,yjl8|.) We are not 
aware of any direct determination yet of the exponent (c in the Edwards- Anderson model. The 
corresponding zero temperature exponent ( has in d = 3 the Migdal-Kadanoff value ( ~ 0.75 



[ [ll| , |l2yiq] ; the difference ( — Cc, just as the difference yc — y, measures the distance to the lower 
critical dimension. The higher moments < AciAej >'^'i and < Ae.^ Ae, >2<;-i decay with powers 



of \i — j\ that involve a hierarchy of exponents C,c{q), see equation (2.38). 



A quantity that should be measurable with relatively good precision is the fourth energy 



cumulant < AE^ >cum (equation (2.41)) in an arbitrary subvolume of linear dimension less 



than the correlation length. As V becomes large, 



<_AE^>^v-^oo \ V^^-^ if 4Cc/d>l. (a) ^^^^ 

^ y est. — (terms nonanalytic in V) if -iQc/d < 1 . (b) 

The MKRG values indicate that in d = 3 the second scenario will be realised, since numerically 
AC,c/d—l = 4-0.57/3 — 1 = —0.24. However, we cannot be sure of the sign of this quantity in d = 3. 
At di, MKRG gives for the same quantity the value of approximately 4-0.75/2.5—1 — 0.2. Hence 
the fourth energy cumulant in d = 3 provides a - possibly rather sensitive - way of determining 

Cc- 
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